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ABSTRACT 

The [Cn] fine-structure line at 158//m is an excellent tracer of the warm diffuse gas in the ISM and the interfaces between molec- 
ular clouds and their surrounding atomic and ionized envelopes. Here we present the initial results from Galactic Observations of 
Terahertz C + (GOTC+), a Herschel Key Project devoted to study the [Cn] fine structure emission in the galactic plane using the 
HIFI instrument. We use the [C n] emission together with observations of CO as a probe to understand the effects of newly-formed 
stars on their interstellar environment and characterize the physical and chemical state of the star-forming gas. We collected data 
along 16 lines-of-sight passing near star forming regions in the inner Galaxy near longitudes 330° and 20°. We identify fifty-eight 
[Cn] components that are associated with high-column density molecular clouds as traced by 13 CO emission. We combine [Cn], 
I2 CO, and 13 CO observations to derive the physical conditions of the [C n]-emitting regions in our sample of high-column density 
clouds based on comparison with results from a grid of Photon Dominated Region (PDR) models. From this unbiased sample, our 
results suggest that most of [Cn] emission originates from clouds with Ht volume densities between 10 35 and 10 55 cirT 3 and weak 
FUV strength (xo = 1-10). We find two regions where our analysis suggests high densities > 10 5 cm" 3 and strong FUV fields 
Cfo = 10 4 - 10 6 ), likely associated with massive star formation. We suggest that [Cn] emission in conjunction with CO isotopes is 
a good tool to differentiate between regions of massive star formation (high densities/strong FUV fields) and regions that are distant 
from massive stars (lower densities/weaker FUV fields) along the line-of-sight. 
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1. Introduction 

The study of the processes governing the formation and destruc- 
tion of molecular clouds is critical for our understanding of how 
galaxies have evolved in our Universe. In terms of column and 
local volume densities only two extreme states of cloud evolu- 
tion have been systematically obse rved: diffuse atomic cl ouds 
traced by the 21 cm line of H i (e.g. iKalberla & Kerp||2009l) and 
dens e molecular cloud s traced by rotational transitions of CO 
(e.g. lDame et al.ll200ll) . We know, however, very little about the 
intermediate phases of cloud evolution and the interface between 
diffuse and dense molecular gas. 

Galactic Observations of Terahertz C + (GOTC+), a 
Herschel Key Project, is devoted to study the [Cn] emission 
in different environments in our Galaxy. The survey will ob- 
serve the [Cn] 158 fim line over a volume weighted sampling 
of 500 lines-of-sight (LOS). Upon completion, it will provide 
a database of a few thousand [C n]-emitting clouds distributed 
over the entire Galactic plane. 

The [Cn] fine structure line at 158 yum is an excellent tracer 
of the interface between diffuse and dense molecular gas. The 
densities and temperatures in this interface allow effective col- 
lisional excitation of this line. The Hi and H2 volume densities 
are a significant fraction of, or comparable to, the critical den- 
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sities for collisional excitationj (3.3xl0 3 and 7.1xl0 3 crrT 3 at 
T = 100K, respectively), the kinetic temperatures are ~100K, 
and the formation of CO is inhibited by limited shielding against 
far-ultraviolet (FUV) photons and therefore most of the gas- 
phase carbon is in C + and some C°. 

Here we present the first results on the molecular cloud- 
atomic cloud interface from the GOTC+ project. During the 
Herschel Priority Science and Performance Verification phase, 
we have colle cted data along 5 LOS s near / = 340° and 9 LOSs 
near / = 20° dVelusamv et al.ll2010h . The focus of this letter is 
to study [Cn] components towards clouds that have sufficient 
column density to have significant 13 CO emission. Such regions 
can be considered as dense Photon-Dominated Regions ( or 
photodissociation regions, or PDRs). PDRs are regions where 
the chemistry and thermal balance is dominated by the effects 
of FUV photons from young stars dHollenbach & Tielensll 19991 
and references therein). These data are therefore important for 
the study of the stellar feedback of newly formed massive stars 
in their progenitor molecular cloud. We combine the [C 11] data 
with observations of 12 CO and 13 CO from the ATNF Mopra 22- 
m telescope to study 58 high-column density PDRs. We use 
the [C n]/ 12 CO and [C n]/ 13 CO integrated intensity ratios to con- 
strain physical conditions of the line-emitting gas comparing 
with a grid of PDR models. 



1 Electrons are a possibly significant collision partner of C + . 
However, the critical electron density for these particles to produce sig- 
nificant [Cn] emission is 9.2 cur 3 at T = 100 K. The density in diffuse 
regions where the abundance relative to H is X(e) as X(C + ) a: 10" 4 is 
modest (~ 10 2 cirT 3 ), while in the denser regions the ionization is sig- 
nificantly lower. In either case, the excitation by electrons is negligible. 
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Fig. 1. [C n], 12 CO, 13 CO, and C 18 spectra towards the line-of- 
sight G337. 826+0.0. The arrows indicate the [Cn] components 
that have 13 CO counterparts. 



The Galactic plane has been studied in [C n] w ith low veloc- 
ity an d spatial res olution observations w ith COBE dBennett et al.l 
11994 and BICE dNakagawaet al.ll998l) . The high angular (12") 
and velocity (0.2kms~') resolution of the Herschel/HIFI obser- 
vations allow us to study for the first time the rich structure 
of molecular clouds along the line-of-sight towards the galactic 
plane. The Kuiper Airborne Observatory allowed the study of a 
han dful of H n regions with velocity resolved fC n] observations 
(e.g. lBoreiko et al.lll988UBoreiko & Betzlll99ll) . However, they 
were limited to massive star-forming regions with dense and hot 
PDRs. The sensitivity of our observations allow us to study for 
the first time the population of PDRs in our galaxy that are ex- 
posed to weaker FUV radiation fields. 



2. Observations 

2.1. Herschel Observations 

We observed the [Cn] 2 P 3/2 -+ 2 Pi/2 line at 1900.5469GHz 

towards 16 LO Ss in the Galactic plane with the HIFI 

dde Graauw et alj|2010h instrume nt on board the Herschel space 
observ atory dPilbratt et al.l 120101) . We refer to IVelusamv et alj 
d2010l) for more details about the [C n] observations. In FigureQ] 
we show sample LOS spectra of [C n], 12 CO, 13 CO, and C ls O. 



2.2. Mopra Observations 

We observed the J = 1 -^ transitions of 12 CO, 13 CO, and C ls O 
toward the observed [C n] LOSs. These observations are part of 
a survey of all GOT C+ positions towards the inner Galaxy be- 
tween / = -175.5°and / = 56.8°conducted at the ATNF Mopra 
Telescope. The angular resolution of these observations is 33". 
Typical system temperatures were 600, 300, and 250 K for 12 CO, 
13 CO, and C 18 0, respectively. To convert from antenna to main- 
beam temperature scale we use a main-beam efficiency of 0.42 
dLadd et al.ll2005h . All lines were observed simultaneously us- 
ing the MOPS spectrometer in zoom mode. The spectra were 
smoothed in velocity to 0.8 km s _1 for 12 CO and 13 CO and to 
1.610ns 1 for C 18 Typical rms noise is 0.6 K for 12 CO and 
0.1 K for both 13 CO and C 18 We checked pointing accuracy 
every 60 minutes using the closest and brightest SiO maser. 



3. [Cn] components associated with molecular 
clouds 

We identify a total of 146 [Cn] velocity components in the ob- 
served LOSs. From this data set we identify components that 
are associated with high-column density molecular gas by look- 
ing for 13 CO counterparts. We identified most of the high- 13 CO 
column density [C n] components by fitting Gaussian functions 
defined by fitting the corresponding 13 CO spectra. The only ex- 
ception was G337. 826+0.0 for which we calculated the inte- 
grated intensity by determining the area within the FHWM of 
the 13 CO emission, as this line-of-sight shows complex veloc- 
ity structure. Based on the 13 CO line parameters we identify 58 
[Cn] components associated with dense molecular gas. All of 
them also show 12 CO emission while 12 show C 18 emission. 
The remaining diffuse atomic and/or diffuse molecular [Cn]- 
e mitting clouds t hat do not have 13 CO co unterparts are discussed 
by lLangeretai] d2010h and I Velusamy et all d2010h . 

In the left panel of Figure |2j we summarize the observed 
characteristics by plotting the [Cn]/ 12 CO and [Cn]/ 13 CO inte- 
grated intensity ratios for the identified components as a function 
of the [C n] integrated intensity. The ratios are calculated from 
integrated intensities in units of Kkms -1 . The mean value and 
standard deviation are 0.29 and 0.6 for the [C n]/ 12 CO integrated 
intensity ratio and 1.75 and 2.54 for [Cn]/ 13 CO The ratios vary 
over 2 orders of magnitude suggesting a wide range of physical 
conditions in our sample. 

We use the [Cn]/ 12 CO and [Cn]/ 13 CO integrated intensity 
ratios to constrain the physical conditions of the line-emitting 
gas. The 12 CO emission, which becomes optically thick quickly 
after a modest fraction of the gas-phase carbon is converted 
to CO, is not very sensitive to the FUV radiation field, as the 
temperature at the C + /C°/CO transition layer is also insensi - 
tive to this quantity dWolfire et al.lll989t iKaufman et all 1 19991) . 
Therefore, the [C n]/ 12 CO ratio is determined by the column den- 
sity of C + and the temperature at the surface of the PDR, which 
are in turn dependent on the FUV radiation field and H2 density. 
The [C n]/ 13 CO ratio is proportional to the ratio between the C + 
and 13 CO column densities. It therefore gives, provided that ex- 
tra constraints on the total column of material are available and 
that there are no significant variations of the FUV field within 
the beam, a constraint on the location of the C + /C°/CO transition 
layer which in turn depends on the strength of the FUV field and 
H2 density. 



4. Comparison with PDR model calculations 

We compare the observed [Cn]/ 12 CO and [Cn]/ 13 CO integrated 
intensity ratios with the results of a PDR model grid in order to 
constrain physical conditions of the [C n]-emitting clouds. 

The model grid was ca l culated usi ng the KOSMA-r PDR 
model dStqrzer et alj 119961; iRollig et all 120061) which is avail- 
able onlinqj- The model provides a self-consistent solution of 
the chemistry and thermal balance of a spherical cloud, with a 
truncated density profile, which is illuminated isotropically by 
a FUV radiation field. The density distribution has the form, 
n{r)-n s {rlr c y is for Q.2r c < r < r c and a constant density of 
n(r) - « s (0.2) in the central region of the cloud (r < Q.2r c ). 
Here r c is the cloud radius and n s is the density at the cloud 
surface. Note that with a power-law index of 1.5, the average 
density of the clump is about twice the density at the cloud sur- 
face. The line intensities are calculated using a non-LTE radia- 
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Fig. 2. (left panel) The observed [C n]/ 12 CO and [C n]/ 13 CO integrated intensity ratios for our sample of [C n] components associated 
with high-column density molecular clouds. The ratios are calculated from integrated intensities in units of K km s~ ' . The horizontal 
error bars are the uncertainties in the determination of the [C n] integrated intensities. The vertical error bars are the uncertainties in 
the line ratios which are derived from error propagation, {central and right panels) Results of the comparison between [C n]/ 12 CO 
and [C n]/ 13 CO ratios for all identified [C n] components and the PDR model grid showing the constrained ranges in FUV radiation 
field (central panel) and H2 volume density (right panel). The arrows indicate lower limits to the H2 volume density. 



tive transfer code by lGierens et al.l dl992l) . Each model is charac- 
terized by the clump mass, the density at the cloud surface, and 
strength of the FUV field. The clump mass ranges from 10~ 2 to 
1OOM , the density at the cloud surface from 10 3 to 10 6 cirT 3 , 
and th e strength of the FUV field from^o = 1 to 10 6 (in units 
of the lDraindll978l fieldl). We do not use Hi and C ls O obser- 
vations to constrain our solutions as model grids involving their 
emission are not available. 

By using a spherically symmetric model we assume that 
the cloud spatial structure can be described by an ensemble of 
clumps with sizes much smaller than the resolution of our obser- 
vations. Additionally, we assume that each clump in this ensem- 
ble has the same mass and density, and that the [C n]/ 12 CO and 
[C n]/ 13 CO line ratios can be estimated using the line ratios of a 
single clump of that mass and density. Therefore, the compari- 
son with the PDR model grid provides the typical incident FUV 
field, mass, and density of the regions that dominate the observed 
line ratios. An even more realistic model considers clumps fol- 
lowing the distributio n of masses and sizes observed in man y 
molecular clouds (e.g. lZielinskv et al.l2000tlCubick et~aT1l2008l) . 

The central and right panels in Figure [2] show a summary of 
the constrained IT volume densities and FUV radiation fields 
for our sample. We consider models with chi-squared^ 2 smaller 
than 1 . \x 2 ■ ■ We find two [C 11] components with high volume 
densities (> 10 5 cm -3 ) and strong FUV fields (between xo = 
10 4 -10 6 ). Both regions are characterized by [C n]/CO integrated 
intensity rati os that are larger than 1 (c.f. Orion has a [C n]/ 12 CO 
ratio of 1.36: ICrawford et al.ll985l) . We show an image and [C n] 
spectrum of one such source in Figure[3] The remaining compo- 
nents have lower volume densities between 10 35 - 10 55 cm -3 . 
Six of those could have a strength of FUV field as high as 100 
while the majority (51 components) have FUV fields between 1 
and 10. For all components the comparison with the PDR model 
grid suggests clump masses that are larger than 1 M . Note that 



3 The average FU V inte n sity of the local ISM is 
2.2xl(T 4 ergcnr 2 s~ 1 sr 1 dDraind [l978h - Note that the Draine 
field is isotropic (i.e. a given point is illuminated from 4n steradians) 
while the surface of the clouds considered here are only illuminated 
from 2n steradians. Therefore the rate of photoreactions at the cloud 
surface are half of what they would be with the Draine field in free 
space. 



due to the limited spatial coverage of the observations presented 
here, the distribution of physical conditions is not smooth. We 
will obtain a much better sampling of the distribution of physi- 
cal conditions in velocity components distributed over the entire 
galactic plane with the completed GOT C+ survey. 

The large number of components with low-FUV field is a re- 
sult of the low observed [C n]/ 12 CO ratios of about 0.1. Such val- 
ues of the [C n]/ 12 CO ratio are expected for^o < 10 3 over a large 
range of H2 volume densities (see e.g. Figure 9 in lKaufman et al.l 
1 19991) . The [C n]/ 13 CO ratio provides an additional constraint on 
the FUV field. The majority of the observed components have 
small ratios that suggest a large column density of I3 CO relative 
to that of C + . This result is suggestive of a C + /C°/CO transition 
occurring close to the surface of the cloud which is a result of ei- 
ther high densities or weak FUV fields. Note, however, that using 
[Cn] and 13 CO to constrain the location of the C + /C°/CO tran- 
sition requires extra constraints on the total column density of 
material thoughout the clump which in turns depends on the as- 
sumed clump surface density and mass. These two quantities are 
not well constrained in the analysis presented here. Additionally, 
it requires that there are no significant variations of the FUV 
field within the beam as shielded clumps might contribute sig- 
nificantly to the 13 CO emission but little to that of [Cn]. The 
H2 volume density for individual velocity components can be 
better determined using observations of the 609pm and 370/im 
transitions of neutral carbon, which have been used to constrain 
the temper ature and density at the C + /C°/CO transition region in 
PDRs (e.g. lStutzkietal.lll997l) . The [C n] to bolometric infrared 
flux is also useful to c onstrain the FUV field dWolfire et al.ll989t 
iKaufman et al.l 19991) but is only useful towards LOSs with a sin- 
gle velocity component. 



5. Discussion 

We have found that most of the [C 11] velocity components con- 
sidered here are associated with regions that are exposed to 
weak FUV radiation field (xo < 10) and therefore are away 
from OB associations. PDRs exposed to weak FUV radiation 
fields ha ve been studied in a few sources using [Ci] emis- 
sion (e.g. iMaezawa efalll 19991; lBenschll2006l; iPineda & Benschl 
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Fig. 3. Example of [C n] emission associated with a massive star 
forming region. The line-of-sight G345. 65+0.0 passes near sev- 
eral bright Hn regions as shown in the Spitzer 8 /mi image. 
In this line-of-sight we find a bright [C n] velocity component 
at -12.9 kms~' which our analysis suggests arises from a re- 
gion with high density (~ 10 57 crcT 3 ) and strong FUV field 

to = io 4 - 5 - io 6 ). 



|2007|) but never observed with velocity-resolved [C n] before 
He rschel/HIFI. 

ICubick et al.l (12008) suggested that most of the [C n] in our 
Galaxy originates from a clumpy medium exposed to a FUV 
field of xo - IO 1 ' 8 , which is larger than the upper limit de- 
termined for the majority of the observed components. Note, 
however, that their model does not consider emission arising 
from diffuse clouds. From our observed LOSs, we find that 
about 56% of the total detected [Cn] emission arises from 
low-column density regions (without sign ificant 13 CO emission; 
lLanger et aLll2Q10tl Velusamy et alJl201(H) while 44% is emitted 
from den se PDRs. Neverthele ss, the moderate FUV field pre- 
dicted by ICubick et al.l (120081) might suggest that the predomi- 
nance of low-FUV radiation field regions observed in our lim- 
ited sample (covering less than 2% of the entire GOT C+ sur- 
vey) might hold for the entire Galaxy. 

We found two regions where our analysis suggests high den- 
sities (> IO 5 cirT 3 ) and strong FUV fields (between^ = IO 4 and 
IO 6 ). These regions are likely associated with massive star for- 
mation. This conclusion is a result of the elevated [C n]/CO ratio 
observed towards these regions. This identification suggests that 
the [C n]/CO ratio is good tracer of the location of massive star 
formation regions in the galaxy. [C n] observations will therefore 
provide an alternative method to determine th e distribution of 
massi ve star forming regions in the galaxy (e.g. Bronfman et al. 
|2000|) . Note that velocity resolved observations are crucial for 
the interpretation of the [Cn]/CO ratio. In our observed LOSs 



we have found velocity components showing [C n] emission but 
no CO as well as components showing CO but no [C n] . Velocity 
unresolved observations would have given a distorted value of 
the [C n]/CO ratio that would result in an incorrect interpretation 
of the physical conditions of the line-emitting gas. 

6. Conclusions 

We have presented velocity-resolved observations of [C n] to- 
wards 16 LOSs located near / = 340° and / = 20° in the Galactic 
plane using the HIFI instrument on board the Herschel space 
observatory. We identified a total of 146 different [C 11] velocity 
components. In this letter we analyzed a sample of 58 compo- 
nents that are associated with high-column density molecular 
gas as traced by 13 CO emission. These components contribute 
44% of the total observed [Cn] emission implying a signifi- 
cantly larger amount of [C n] emission originating in the diffuse 
ISM than from star forming environments. We have compared 
the [Cn]/ 12 CO and [Cn]/ 13 CO integrated intensity ratios with a 
PDR model grid to constrain the strength of the FUV field and 
the H2 volume density in this sample. We find two clouds for 
which our analysis suggests high densities (> 10 5 cm~ 3 ) and 
strong FUV fields (xo — IO 4 - IO 6 ), likely associated with mas- 
sive star formation. The majority of the observed components, 
however, have modest densities (10 35 - 10 5 5 crrT 3 ) and weaker 
FUV fields (xo - 1-10). Although the population of clouds with 
these conditions is likely where most of the [C 11] emission origi- 
nates in our Galaxy, their properties are largely unexplored. The 
GOT C+ survey will provide a few thousand clouds distributed 
in the Galactic plane and therefore will allow us to characterize 
this population of intermediate clouds. 
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